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Abstract — In this paper, a mathematical model that describes 
the interaction of a pulsed mode ultrasound waves and a 
multilayered biological structure is developed. As a special case, 
the model is applied for quantitative monitoring of fractured 
bone healing. Four different frequencies are considered. It is 
found that a transducer operating at 1 MHz has the highest 
sensitivity to the different phases of callus hardening. On the 
other hand, a 3 MHz transducer yields the highest sensitivity to 
the thickness of the callus layer. Thus, using both frequencies 
might prove to be useful for full quantitative monitoring of the 
bone healing process. 

Index Terms — Ultrasound, Pulsed mode. Bone fracture. Tibia. 

I. INTRODUCTION 

When a bone is fractured, treatment involves reduction or 
realignment of the bone, immobilization and restoring 
function through rehabilitation. The healing process starts 
when new bone (callus layer) forms and fills in the fractured 
area to restore the original continuity and solidity. It generally 
goes through the same series of stages for all fractures: 
inflammation, soft tissue formation, hard tissue formation, 
new bone and remodeling [2]. As healing progresses, a soft 
callus layer forms and begins to harden with time. At a later 
time stage, the callus layer begins to decay. Consequently, by 
monitoring the hardening and decay of the callus layer, bone 
healing stage may be quantitatively estimated. Healing is 
considered to be complete when the bone is almost identical to 
its original shape before injury and has regained its normal 
stiffness and strength. A bone fracture may be diagnosed by 
many diagnostic imaging tests such as X - Rays, Computed 
Tomography (CT), Magnetic Resonance Imaging (MRI) or 
Vibratory devices [3-5]. X - Ray radiation is hazardous and is 
unsafe for repetitive tests. Moreover, fractures are typically 
visible on primary radiographs. However, the healing process 
during the early soft tissue formation stage is difficult to 
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visualize on radiographs. Furthermore, X-rays can not be used 
for pregnant women or patients who had a barium contrast 
media or any medication containing bismuth. MRI tests are 
long and of high cost and can not be repetitive. The MRI 
machine itself is prohibitively expensive for small hospitals 
and therefore, it is not available everywhere. Vibratory 
devices are mechanical sine wave vibrators that are applied to 
the fracture site [5]. An electromagnetic shaker is pushed 
against the skin near the distal end of the bone. A miniature 
accelerometer captures the response of the bone-prosthesis 
system to the applied vibrations. If there is no fracture, the 
bone system behaves linearly and only the excitation 
frequency will be detected. However, if there is a fracture, the 
system will behave in a non-linear mode and harmonics will 
be detected by the accelerometer. Experiments have shown 
that an excitation frequency between 100 Hz to 200 Hz is well 
suited to detect fractures. Unfortunately however, vibratory 
devices suffer from the increased signal damping and locating 
the resonance frequency especially with obese subjects. 
Recently, the potential of utilizing Ultrasound was 
demonstrated as a tool to monitor the healing process of the 
Tibia using continuous mode ultrasound [6]. In this paper, 
pulsed mode ultrasound is assumed. A mathematical model 
that is suited for practical pulsed mode applications is 
developed. The model is applied to the problem of monitoring 
of the tibia bone healing process. The paper is organized as 
follows. Section II presents the theoretical derivation of the 
mathematical model. Results are presented in section III. 
Finally, conclusions are drawn in section IV. 

II. THEORY 

In this model it is assumed that an ultrasound transducer is 
radiating into an N-layer structure, terminated by an infinite 
half space. Fig. 1 depicts a 3 -layer structure terminated by an 
infinite half space of any material with acoustic impedance Z4. 
Each layer is characterized by four parameters: 

Thickness of the i^^ layer in (m), df 
Acoustic impedance of the i layer in (kgm" s" ), Z/ 
Velocity of the wave in the i^^ layer in (ms"^), Vt 
Attenuation factor of the i^^ layer in (dBm"^), a. 

According to the reflection laws of ultrasound waves, waves 
experience reflection whenever they travel from one medium 
to another. The ratio of the pressure of the reflected wave to 



PROC. CAIRO INTERNATIONAL BIOMEDICAL ENGINEERING CONFERENCE 2006 



1 



the pressure of the incident wave is called the reflection 
coefficient (F/) and it is given as 



r/ = 



z+l / 



In addition to the main first reflections at each interface, there 
exists additional multiple reflections within each layer. The 
first and second multiple reflections and their corresponding 
times may be obtained in a similar fashion to the main 
reflections. The first set of multiple reflections appearing at 
the surface of the transducer occur at: 



Where /] is the reflection coefficient from the interface 
between the {^ layer and the (i+1)^^ layer. 
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The second set of multiple reflections appearing at the surface 



Fig. 1: A 3 -layer structure terminated by an infinite half space 

of a material. Each layer / has thickness di, acoustic 
impedance Z/, propagation velocity F/ and attenuation factor 
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For pulsed mode ultrasound, it is assumed that the incident 
wave is a delta function at time t=0. Fhe wave travels a 
distance di before reaching the first discontinuity. At the 
interface, the wave will split into two components: one will 

propagate to the next layer with amplitude (1 - Fi ) and the 

other will be reflected and will travel a distance of di to the 
transducer with amplitude F^ . Fherefore, the first reflection 

which is equal to F. , will be recorded at the transducer end at 

Id^ 
the elapse of time . Fhe amplitude of the second 



reflection is equal to (1 - Fi ^Yn. (1 + Fi ) which will be 
recorded at the transducer end at the elapse of 



time + . Finally, the amplitude of the third reflection 

is equal to (l-F^ )(1-F2)F3(1 + F2)(l + Fi ) and will be 
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b) time t = + + with amplitude 

F3 V^ V, 

(1-Fj)(l-F2)F3(-F2)F3(-F2)F3(1 + F2)(l + Fj) 



Although multiple reflections will continue to take place, 
usually the first two sets of multiple reflections are of 
significant amplitude relative to the first reflection. Therefore 
in this work the first two sets of multiple reflections will be 
considered. 

A unit impulse was assumed to be the input in the derivations 
outlined earlier. However, to simulate a practical setup, a 
dying oscillatory pulse with a finite duration should be 
considered. Hence, the output is obtained by convolving the 
impulse response with the oscillatory input pulse /fifj. 
The main and multiple reflections resulting from the structure 
mentioned earlier show a recursive pattern. Therefore, this is 
exploited to get a general formula for an N-layer structure. 
The main reflection in each layer Pi(i,t) and the multiple 
reflections ^2^?,^^^ are: 
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Where the attenuation coefficients are: 
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III. RESULTS 



^. Model Validation 



To vahdate the proposed model, an Aluminum plate 
immersed in water was used in a laboratory setup. An 
ultrasound pulser receiver made by OPTEL was used in 
conducting an A-scan [7]. The structure parameters are as 
follows: 



St 



-2 -1 



1 Layer: Water, Z = 1.5x10 kgm s , d= 1 
cm, V= l4S3ms~^ , a= 0.002 dBcm~^MHz~^ 



,nd 



-2 -1 



• 2 Layer: Aluminum, Z = 17x10 kgm s 

1 cm, V= 6320w^"^ , a= 0.0182 dBcm~^MHz~^ 

• 3^^ Layer: Infinite half space of water. 

Fig. 2(a) shows the measured incident and reflected voltage 
as a function of time. The measured signal includes the 
main and multiple reflections. The first three prominent 
bursts originate from the transducer, while the rest of the 
reflections correspond to the main and multiple reflections 
of the structure. The used measurement system included a 
time-variable gain amplifier to compensate for the 
ultrasound wave attenuation as a function of time [7]. 
The peak of the second burst coincides with the time at 
which the incident wave leaves the surface of the 
transducer. Consequently, the first main reflection due to 
the water aluminum interface appears after 13 jis (it appears 
at t ~ 17 |is). The second main reflection due to the 
aluminum water interface appears approximately at 20 |is. 
The obtained results indicate that the measurement system 
is properly calibrated. 

The developed model was tested using the same multi 
layered structure described above. Fig. 2(b) shows the 
reflected pressure wave as a function of time. These waves 
include the main and multiple reflections. The first three 
bursts are assumed to originate from the transducer with 
time referenced to the middle of the second burst (to mimic 
the experimental results) while the rest of the reflections 
correspond to the main and multiple reflections of the 
structure. Consequently, the first main reflection due to the 
water aluminum interface appears after 13 jis (it appears at t 
-17 jis). The second main reflection due to the aluminum 



water interface appears approximately at 20 jis. The 
obtained results reflect good agreement between the model 
and real measurements. The simulated results show the 
effect of ultrasound attenuation as a function of time. The 
time-variable gain amplifier was not included in the 
simulation. 
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Fig. 2: Pulsed ultrasound A-Scan of a 3-layer structure 
composed of water-aluminum-water, 
a) Experimental b) Simulation. 
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B. Biological Structure 

The model described above was used to model the 
ultrasound interaction with biological structure composed 
of skin, fat, callus, and infinite half space of bone as shown 
in Fig. 1. The structure parameters are as follows: 
• Layer 1: Delay line of water, 

thickness (di)= 6 mm. 



Z = 1.5x10 kgm s 



Velocity (Vi)= \A^3ms , attenuation factor 

(ai)=0.002dBcm~^MHz~^ . 

Layer 2: Skin, Z = 1.69x10 kgm s , 

d2=2.5mm, V2= 1531ms , a2= 

3.5dBcm MHz 

Layer 3: this layer is a fat layer with a thickness 
(ds) of 10 mm except for a small circular area (1 
cm^) where the fat is replaced by callus. The fat 

layer has Z = 1.38x10 kgm s 



Velocity= 

1450ms , and attenuation= l.SdBcm MHz 
The callus layer has a thickness of 10 mm with 
properties similar to fat at the early stages of callus 
formation. The callus acoustic properties change 
from those of soft callus (similar to fat) to hard 
callus (similar to bone). In other words, its 
acoustic impedance, velocity of the wave and 
attenuation will increase with time (healing). After 
the callus hardens, it will begin to decay from 10 
mm to mm indicating complete healing. Both the 
callus thickness and its properties depend on the 
healing stage. Consequently, monitoring these 
quantities will reveal information about the 
healing stage. 

/^ 1 

• Layer 4: Bone, Z = 7.8x10 kgm s , d4 is an 
infinite half- space (usually 3 mm thick), ¥4= 

3500ms~ , a4= \3dBcm~ MHz~ . 

Four different frequencies were considered in the 

simulation: 1 MHz, 3 MHz, 5 MHz and 10 MHz. Fig. 3 

shows the reflection coefficient variation as a function of 

callus acoustic impedance changes for different callus 

thicknesses. The simulation results show the remodeling 

stage of the bone healing process. In Fig. 3a, for example, 

each curve corresponds to a certain callus thickness starting 

with soft callus of acoustic impedance similar to fat and 

ending with hard callus, with acoustic impedance similar to 

bone. A reflection coefficient of 0.3 represents acoustic 

^x ^kg^_2^_^ 

WJfi§^^^^ocfi^i^n^ithl^M^%ict^te^ l^ween to 10 

mm. This uncertainty may be resolved using the A-scan 

which reveals the thicknesses of all layers. Fig. 3b, c, and d 

show similar results using sensors operating at different 

frequencies. 

Our criterion to compare the performance of the different 

transducers is the difference in the magnitude of the 



reflection coefficient, which reflects the sensitivity of the 
transducer to different stages of bone healing. As can be 
seen, the most sensitive frequency to acoustic impedance, 
and hence the properties of the callus is 1 MHz since it 
produces a maximum reflection coefficient difference of 
0.55 for callus thickness of 10mm. This difference 
corresponds to the change from soft callus to hard callus. 
However, if we target the thickness of the callus, using a 3 
MHz sensor was found to be optimum since it produces the 
maximum difference in the magnitude of the reflection 
coefficient due to different callus layer thicknesses. This 
suggests that multi frequency measurements, using 1 MHZ 
and 3 MHz, are needed to reliably determine the status of 
the callus and its thickness. 
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Fig. 3: Magnitude of the reflection coefficient as a function 
of the acoustic impedance of the callus layer at 
different thicknesses operating at a frequency of 
a) 1 MHz. b) 3 MHz. c) 5 MHz. d) 10 MHz. 

Fig. 4 shows the position of the reflected wave as a 
function of the thickness of the callus layer. We remark that 
this measurement parameter is independent of the 
frequency and the status of the callus, as can be interrupted 
from the mathematical model. The position of the reflected 
wave changes linearly with the thickness of the callus. This 
is extremely important from the practical point of view. It 
means that quantifying the thickness of the callus is 
possible due to the linear behavior. Additionally, since 
different people have different thicknesses of their limb 
tissues, a shift in time will occur without changing the 
overall pattern observed in Fig. 4. 
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Fig. 4: Position of the reflected wave as a function of the 

callus thickness. 

IV. CONCLUSION 

A mathematical model that describes the interaction of 
ultrasound waves with a multilayer structure was presented. 
Preliminary experimental results agree with the simulation 
results. As a special case, the model was applied to 
monitor bone healing process. A 3 -layer structure backed 
by an infinite half-space was used. Transducers operating at 
frequencies of 1, 3, 5 and 10 MHz were assumed. Results 
show that transducers operating at frequencies of 1 MHz 
and 3 MHz can reliably be used to monitor the bone healing 
process since they provide highest sensitivity to callus state 
and thickness. 
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